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Abstract 



The pseudo-SU(3) model is extended to explicitly include the spin and proton- 
neutron degrees of freedom. A general formalism for evaluating matrix elements of 
one-body and two-body tensor operators within this framework is presented. The 
pairing interaction, which couples different irreducible representations of SU(3), is 
expressed in terms of pseudo-space tensors and a general result is given for calcu- 
Tij- . lating its matrix elements. The importance of pairing correlations in pseudo-SU(3) 

^1* I model calculations is demonstrated by examining the dependence of wavefunctions, 

low- energy collective excitation spectra, and moments of inertia on the strength of 
the pairing interaction. 

'b 

• • : 1 Introduction 

> ■ 

•i-H . 

^ ' Over the past few years there have been a number of discoveries in nuclear structure 

j^ ■ physics (superdeformed shapes, identical bands, scissors mode, etc.) which have raised a 

number of interesting questions, many of which remain unanswered. In addition, one can 
anticipate the identification of additional new phenomena in the near future due to ongo- 
ing improvements and innovations in experimental techniques {Air detectors, radioactive 
beams, etc.). These experimental developments seem not to be matched by a comparable 
record of theoretical achievements focused on the questions that are being raised. 

However, it is commonly accepted that the nuclear shell model should be able to 
address these issues and provide answers to these questions. The problem is that most 
shell-model theories are limited by the large dimensionalities of the required model spaces. 
For example, existing algorithms and available computing equipment still only allow for 
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reasonably complete standard shell- model calculations in the mass range A < 28. And 
even the best of these calculations do not include a sufficient number of multi-shell config- 
urations to adequately address the issues raised by experiments that probe details of the 
momentum and current distributions in nuclei. Furthermore, most of the new discoveries 
are for much heavier nuclei, typically in the mass range A > 100. 

When faced with a situation like the one just described, it is essential to take full 
advantage of symmetries, those which are only approximately fulfilled as well as those 
that are exact. The selection rules that are associated with such symmetries generate, 
respectively, weakly coupled and disconnected subspaces of the full space and this in turn 
allows for a significant (orders of magnitude) reduction in the dimensionality of the model 
space. Since matrix elements of the system's Hamiltonian between these subspaces only 
vanish if the symmetry is exact, the question that must be addressed is whether or not 
using the identified approximate symmetry is physically justified. 

Pseudo-spin |l|, |[ is such a symmetry and it serves as the conceptual foundation for 
the subject of this contribution, namely, the pseudo-SU(3) model [Q, and its extension, 
the pseudo- symplectic model [^. The success of the pseudo-spin concept relies on the 
fact that the harmonic oscillator structure of the usual shell-model Hamiltonian, which is 
destroyed by the spin-orbit interaction, is restored under a transformation to the pseudo- 
space/spin picture. The idea is to express all normal-space shell-model quantities in terms 
of their pseudo-space equivalents by applying a normal <-^ pseudo transformation to both 
operators and wavefunctions [§. Since for practical purposes the harmonic oscillator 
structure of the Hamiltonian is restored in the pseudo-space scheme, this transformation 
makes it possible to invoke powerful group theoretical truncation methods so shell-model 
calculations can be carried out even for heavy nuclei with A > 100. 

There are still questions being raised about the physical implications of the approxima- 
tion one is making by assuming that the pseudo-spin symmetry is exact. However, there 
have been a number of publications p, ^, ^, ^ which demonstrate in both single-particle 
and many- particle systems that even if it is not an exact symmetry, for practical appli- 
cations pseudo-spin is a sufficiently good symmetry to be of major physical importance. 
The reader is referred to these publications for assurance regarding the validity of the 
pseudo-spin concept. Here we focus only on an algebraic many-particle shell-model the- 
ory that takes full advantage of pseudo-spin symmetry, namely, the pseudo-SU(3) model. 
The usefulness and predictive power of this model have been demonstrated in many ap- 
plication ranging from the calculation of collective excitation spectra and a scissor's 



mode study [|TU[ to fundamental problems likes the half-lives of /3/3-decay modes |TT . 

In this contribution we propose extensions of the pseudo-SU(3) model in order to 
widen its range of applicability and in so doing increase the capabilities of the model for 
describing and predicting new and interesting experimental phenomena: 



The ansatz of reference ||I0[ is generalized by explicitly including spin degrees of 
freedom in a full proton- neutron formulation of the theory. This generalization 
makes it possible to describe properties of even-odd, odd-even, and odd-odd nuclei 
in addition to those of even-even ones. Within this extended approach, tensor 
operators are defined and expressions for calculating their matrix elements are given. 



In spite of many successful applications of the pseudo-SU(3) model, the Hamiltonian 
that has been used in most applications is highly schematic in nature and needs to 
be extended. The usual structure taken for the pseudo-SU(3) model Hamiltonian is 

H = H,-^Q^-Q- + i?,oii (1) 

where Hq is the harmonic oscillator Hamiltonian, Q°^ denotes the algebraic quad- 
rupole operator (see ref. 0]), and -R^qH stands for a residual interaction that is 
introduced to describe properties of low-energy collective excitations. A comparison 
with other shell- model theories indicates the importance of introducing a short range 
interaction, and therefore it seems reasonable to add a pairing interaction to terms 
already included in Eq. (|l]) [0, [l^ . 



The reason the pairing interaction has not been introduced previously is a technical 
matter: Basis states of the pseudo-SU(3) model are organized according to irreducible 
representations (irreps) of SU(3) which are labeled by (A, /x). In the past it was not possi- 
ble to calculate matrix elements between different SU(3) irreps (except for a few relatively 
simple cases), but recently a general code was released that removes this limitation fl^ . 



The point is that the terms that are part of the traditional Hamiltonian of Eq.(|^) do 
not couple different SU(3) irreps while the pairing interaction has non-vanishing matrix 
elements between states of different (A, /x). It follows from this that the calculation is far 
more interesting and challenging when the pairing is present [|1^ . 

The paper is organized as follows: In the next section the structure of basis functions 
and operators of the generalized pseudo-SU(3) model are defined and general expressions 
for evaluating one-body and two-body matrix elements are determined. In Section 3, the 
pairing interaction is expanded in terms of SU(3) tensor operators and expressions for 
its matrix elements are given. The importance of pairing for the pseudo-SU(3) model 
is illustrated in Section 4 by investigating its influence on wavefunctions and excitation 
spectra. A summary is presented in Section 5. 

2 Wavefunctions, operators, and matrix elements 

2 . 1 Wavef unct ions 

Wavefunctions of the pseudo-SU(3) model which take the spin and proton-neutron degrees 
of freedom explicitly into account are labeled by the eigenvalues of the Casimir invariants 
of the group chains depicted in Fig. |I]. 

Figure ^ 

Since the pseudo-SU(3) model is a many-particle theory, it is to be understood that its 



wavef unctions are subject to the Pauli Exclusion Principle, that is, they must be anti- 
symmetric under the interchange of identical particles. Because a particle-permutation 
<-H> unitary-symmetry complementarity applies, this requirement can be met by restricting 
the basis states to the antisymmetric irreps of U{2Qf^), where a = it for protons and a = u 



for neutrons [|T^ and where fio- = {rjcr + 1)(^ct + 2)/2 is the number of normal parity levels 
in the rj^r shell. 

These anti-symmetric irreps of f/(2fio-) are characterized by l^"^ where Nfj is the 
number of normal parity particles. The N^r in this result can be determined by a simple 
procedure with the help of a Nilsson diagram (see, for example, [|l3): At the experimen- 
tally determined P2 value, the proton (neutron) Nilsson levels (unique parity as well as 
normal parity) are filled pairwise (to accomodate the spin degree of freedom) from the 
bottom of the well up until the total number of protons (neutrons) equals the number in 
the nucleus under consideration. If shifts in j32 by 10 to 20 percent of the original value 
do not influence the distribution of particles in the unique and normal parity orbitals, one 
can assume that the Ncr determined in this way is more or less unique for that particu- 
lar nucleus. (A Nilsson-like single-particle code which improves upon and automates this 
procedure is currently being developed. Specifically, for any given nucleus it calculates N^r 
as a function of the deformation /32 together with statistical measures that characterize 
the sensitivity of the A^^- to changes in /?2.) 

The irreps of f/(f2o-) and U{2)^°''> are denoted by [/o-] and [f^], respectively. The overall 
antisymmetry in U{2Qfj) and the particle-permutation ^-* unitary-symmetry complemen- 
tarity means that [fa] is related to [fa] through row <-> column exchange ||T^. Since the 
overall attractive nature of the nucleon- nucleon interaction insures that spatially sym- 
metric irreps of the U{Qa) he energetically lowest, in this work only the most symmetric 
irreps of (/{Q-jr) and U^Q^) are taken into account. The physical consequence of this re- 
striction is that no couplings to spin-flip excitations are considered. This approximation 
can be justified by noting that spin- flip modes normally lie significantly higher in energy 
than the low-energy collective excitations which are of primary interest here. 

The reduction U{fia) ^ SU{3) (see Fig. |I|) determines the (Ao-,/io-) irreps that are 
contained in a given [fa]. Multiple occurences of the same (Ao-,/io-) in a fixed [fa] are 
distinguished by a running integer index Oa- The SU(3) labels of the total wavefunction 
are calculated by taking all possible products {{Xt,,^^^) x {Xu,fJ'u)} —^ (-^jA*) into account 
and numbering multiplicities in this product by a running integer index p. The total 
orbital angular momentum L is determined via the chain SU{3) D 5*0(3) where k labels 
its multiplicity. Likewise, the total spin of the protons, S-,,, and neutrons, S^, are coupled 
to the total spin of the nucleus using the usual rules for coupling angular momentum. And 
finally, the total angular momentum J results from the coupling of total orbital angular 
momentum L with the total spin S. 

These group structures, reductions, and couplings are an integral part of the pseudo- 
SU(3) model and all of this information enters the corresponding wave functions: 

I {A'7r[/7r]«^ (Ayr, ^tt) , A^i.[/i/]a,. (A^,, flu) } P (A, /i)fi:L {S^, Su} S ; JM) 



= X!( (^'T, fJ'Tr) KnL^ML^ ] {X^, /i^) K^L^Ml, \ (A, /i) kLMl )p 

{-} 

X {S^Ms^, S,Ms^ I 5M5) X {LMl, SMs \ JM) 

{-} 

X {S^Ms^, S^Ms^ I SMs) X (LM^, SM^ | JM) 

(2) 

with the abbreviation {-} = {Ms^, Ms^, Ms, Ml^, Ml,, Ml, k^, k^, L^, L^} . In this 
result the {LiMl^, L2ML2 \ LMl) factors are SU(2) Clebsch-Gordan coefficients and the 
SU(3) couphng coefficients (. . . ; . . . | . . .) and (. . . ; . . . || . . .) are defined in |l^. The no- 
tation of Eq.(^ indicates the following couplings: 

• SU(2): [L^ X L^]^*^^ [S^ X 5^]^^'^^ and [L x 5]^^^ where [A x Bf^^^ denotes the 
usual angular momentum coupling [^, pO|| ; 

• SU(3): {(Att, fijr) X (X^, fi,^)} — > {p{X,fi)KLM} where the multiplicity index p 
numbers the multiple occurences of (A,/i) in the product {{Xt,, p^) x (A^,, Pu)}. 

2.2 Operators and their matrix elements 

As explained in the introduction, accepting the fact that SU(3) is a good pseudo-space 
symmetry is a fundamental assumption of the pseudo-SU(3) model. The validity of this 
approximation, which applies to nuclei with A > 100, is based on: 

• the goodness of pseudo-spin symmetry, which is suggested by the Nilsson model and 
confirmed by other more sophisticated theoretical analyses p, H, 0; 

• the observation that most heavy nuclei are deformed, which is important because 
the goodness of the SU(3) symmetry increases with increasing deformation. 

And as shown in Section 2, there is a group-subgroup chain for labeling basis states that 
has SU(3) as an integral part of its structure, and which therefore takes full advantage of 
the pseudo-SU(3) model symmetry. To guarantee wavefunction-operator compatibility, 
it is also necessary to express operators associated with physical observables in terms 
of SU(3) tensors. The procedure for doing this can be found in reference |]10|- In this 
contribution, the spin and proton-neutron degrees of freedom are explicitly incorporated 
into that formalism. For the sake of completeness, it is useful to repeat some of the 



material given in reference |]TOl as it forms a backdrop for the introduction of an extended 



notation that allows pairing correlations to be incorporated into the pseudo-SU(3) model. 
The formalism requires one to know the relevant group's coupling coefficients and the 
tensorial properties of the single-particle creation and annihilation operators, which for a 
fermion system obey the usual anticommutation relations: 

{a+, ak}+ = Sik, {af, a^}+ = {aj, ak}+ = , (3) 

where i (and k) denote the full set of single- particle quantum numbers \ (r/, 0) Im ^nis >. 
Here {ri,0) are the single-particle SU(3) labels (A, /i) and /, |, m and m^ stand, respec- 
tively, for the single- particle orbital angular momentum, spin, and their z-projections. 
The effect of a creation operator on the single-particle vacuum is given by 

^lo)im^^mJ^) = \(v,0)lmlm,). (4) 

While the creation operator at n.^, i is a proper SU(3) tensor, the annihilation operator 
is not and must be modified appropriately to achieve this status 0, [10|, ^ : 



''{0,rj)lm. ^ms = y~^) ^{vfi)l-m ^-nis ■ V^J 

With these elementary definitions in place, it is possible to construct three distinguish- 
able tensor operators as products of two creation and/or annihilation operators: 

• A one-body (unit tensor) operator, 

(A, fj.) kLMl SMs 



(1,1) ^{X,fi)KLMi,SMs ^ / + ~ 



= E((^''0)1^'(0'^)1^II(^'/")'^^) X {l'm',lm\LML) 
{-} 

X ( ^<, ^rris I SMs) a+,_^^^,^, i^, d^o,rj) Imlm. , (6) 

with the abbreviation { — } = {/, m, m^, /', m', m'^} and where the left superscript 
denotes respectively the number of creation and annihilation operators in the prod- 
uct. Note that the multiplicity index p connected with the generic SU(3) product 
{(Ai, /ii) X (A2, /i2)} —^ P (A, /i) is usually dropped when the coupling is mutliplicity 
free, that is, whenever pmax = 1- To summarize, the couplings in Eq.(|^) are: 

- SU(2): [/'x/]^^^ [|xi]^^'^s. 

- SU(3): {(r/', 0)x(0, r7)}^{pmax = l,(A,/i)} . 

An operator that creates a pair of particles, 

r ^ (A, m) kLMl SMs 



where the coefficients that enter to effect the SU(2) and SU(3) couphngs in this case 
are similar to those that appear in Eq.(^. 

• A pair annihilation operator, 

(0,2) TT (m, A) ^LMl SMs - (~ ^ ~ \ ('^' ^) '^^^^^ ^^S 

•^(0,r,')(0,r?)ii - r(0,r,')| ^>,»?)i/ 

In a completely analogous manner, it is possible to define two-body SU(3) unit tensors 
as coupled products of the (^'^^^ and (°'2) J^ operators: 

(2,2) ^pi^,tJ.)KLML SMs ^ /{2,0)t-(Ai'/^i)'S'i ^ (0,2) ^-(^2, A2) S2 V^ ' ^' '^ '^ ^ f^\ 

•^(Ai,/.i)(/.2,A2)5iS2 - \ •^(r,i,0)(7?i,0)ii ^ '^(o,,,^)(o,r,2)|i J " ^^^ 

Note, that the indices r][, r]i, Vj^-, 772 are suppressed in ^^'^^JF. This simplifies the nota- 
tion and can be done because these quantum numbers enter into a pseudo-SU(3) model 
calculation as fixed parameters. In fact, the constraints ti[ = rji, r]2 = r]2, and rji = 772 
apply since inter-shell couplings are not part of a Ohu theory. These restrictions are lifted 
in a pseudo-symplectic extension of the pseudo-SU(3) model Q, and therefore for such 
applications the notation must be expanded appropriately.^ 

The definitions of *^^'^)jF and "^^'^^JF as one-body and two-body operators, respectively, 
are applicable to identical particle systems. The next logical step is the introduction of 
proton- neutron SU(3)-coupled operators. For example, the simplest multipole-multipole 
proton- neutron tensor operator can be defined as follows: 



•^(A.„,M.o)(A.„,M.o)-S.„S.„ = \*^«„,0)(0,7?.Jii «,0)(0,r,.jii 

= Y. {LoMl^.SoMsAJoMo) 

MLoMSo 

"^ 1-^(,,;^,0)(0,7?.jii "^ •^«,0)(0,ry.jii 



po ('^05 1^0) ^oJ-'oOo'i Jo-^io 



p{Xo,flo) t^oLoML^SoMso 



(8) 



The subscript "o" in this expression refers to "operator" and is introduced in view of 
the forthcoming calculation of matrix elements in the pseudo-SU(3) model. The left 
superscript ttz/ on the '""J^ operator, which is a shorthand notation for the complete label 

^A complete theory would also include a general discussion of non-particle number conserving as 
well as particle number conserving operators in both normal-ordered (creation operators to the left of 
annihilation operators) and multipole (for example, products of one-body tensors) forms. Indeed, it should 
be clear that an n-body operator (("'"' J^) can always be expanded in multipole form (and vice- versa, of 
course) by using the fundamental commutation relations, Eq.(||), together with recoupling formulae to 
appropriately rearrange (group/order) the creation and annihilation operators (see |22| ) . Non-particle 
number conserving operators ((™'")^ with m j^ n), represent pickup (stripping) type phenomena when 
m < n (m > n). Here it is enough to limit the development to the forms already introduced as these 
suffice for our Hamiltonian, inclusive of the pairing interaction. 



{(mTT, nT,){my, n^)} = {(1, 1)(1, 1)}, denotes the fact that this operator acts simultaneously 
in the proton and neutron spaces. The left superscript on each of the JF factors in Eq. (^ 
is suppressed because it is (1,1) in both cases, that is, J-' = ^^'^^J^. 
In general a more explicit and complete notation is required: 



7rZ> '^Po\^0; f-to) l^oLiobo'i Jo^'io 



-.)S. 






U 'T'Pi^o i^i^o 5 P't'o 7 '-'f o 



^ } 



Po (AoiMo) KoLoSo; JoMo 



(9) 



where {tt} stands for { [a^,(A^^, /x^J, S'^J x [^^^(/i^a, A^J, S^^] -^ p^„ (A^„, /i^J, S^^} and 
similarly for {z/}, and the a^r and a,^ labels represent additional quantum numbers that 
are required for a unique identification of the '^T and ^JF factors in the product. In 
this case the left superscripts tt, i/, and vrz/ stand for the sets (771^,71,^), {711^,11^), and 
{(m7r,'n.,r)("^iv, '"-£/)}, respectively. The matrix elements of such operators can be calcu- 
lated if the matrix elements of the factors are known: 



{{KIKWAK, f^'.),K[f>'AK, f^'.)}p'{y, f^')'^'L' {S'^,si} s'-, f 



\7VU ^Po (Ao, po) KoLioOo'i Jo I 



{A^Trl/TrlaTT (A,r, AItt) , Ny[fy]ay {K, fi^) } p (A, /i) kL {S^, S^} S ; J) 



x{ S So S' > xl S^ Sy^ SI }>^((A, /i)fi;L,(Ao, /io)fi;o^o|| (A,' /i')«^'^')p 

(7 fj Q U \ I O Oq tJ 

{(A^, /i^) (A^„, /i^j (a;, ^;) p^ ' 

(A^, /i^) {K^, py^) (A'^, ^;,) p, 

(A, /i) (Ao, /io) (A', /i') p 
P Po P' 

^1 v'^'^i ' Pi^\ )^i^2 \^^2 ' ^^2 •' '-''^1 ^^2 



X (iv:[/^]«t(A:., z^:.)^: m ^^::ft:;:::Jtv..,. .. 5.. 111 a^.[/.]«.(a., ^.)^.)p. 



(10) 

where x{- • •} denotes a unitary 9-j or Jahn- Hope symbol |^ ^ and {. . .} its SU(3) 
extension p^. The triple-barred matrix elements that enter in Eq. (|10|) , (. . . ||| . . . ||| . . .), 
are reduced with respect to both SU(3) and S0(3) (see Appendix |B|), and can be evalu- 



ated for operators of physical interest with the code introduced in reference [14|. When 



this result is used to evaluate matrix elements of either a pure proton {{m,^, riy) = (0, 0)} 
or a pure neutron {{m.„,nT,) = (0,0)} operator, it simplifies just as for the corresponding 
SU(2) case; namely, the occurrence of a zero-body factor with its null SU(3) irrep char- 
acter means that the 9-(A,/x) coefficient can be reduced to a 6-(A,/i) coefficient and the 
corresponding triple-barred reduced matrix element goes to unity. The pairing interaction 
is a sum of two such operators: Hp -^ {(2,2){o,o) jr _^ (o,o)(2,2) jr|_ 

3 The pairing interaction 

As stated in the introduction, the inclusion of a short- range interaction of the pairing 
type in pseudo-SU(3) calculations has not been tested and seems from other independent 
evidence to be important. But before exploring the influence pairing has on collective 
nuclear properties, it is necessary to introduce some additional formalism. In second 
quantization the pairing Hamiltonian [I3| is deflned as: 



G 

Hp = — Ty^a~^,i- a^--r- a_ T -^ _ a„ , 1 • ^, , (11) 

^ 4 r1 '''I'l^Ji™! Vih^Jimi '72'25J2'"2 mh^J2m2y V / 

where {-} = {771 , ^i , ji ,"^1 ,^2 , ^2 , J2 ,"^2} with (m^ = -ji,-ji + l, . . . ,ji) and (i = 1,2). 
The bar over quantum numbers in Eq. ([Tl|) denotes time-reversal and G is the pairing 
strength, which is normally taken to be somewhat different for protons and neutrons. 

Obviously Hp, as given in Eg . (|TT1) , is deflned in terms of normal-space quantum num- 
bers. However, since all calculations in the pseudo-SU(3) model are performed in the 
pseudo-space, it is necessary to transform Hp into its pseudo- space representation. An 



application of the normal — > pseudo mapping is described in reference |T0| for the gen- 



eral case, and in reference |^5|] for the special case of pairing. The following is the flnal 



identical-particle result, expressed in terms of pseudo-space quantum numbers [0: 

Hp = -f E(-l)'^^'' ((2^1 + 1)(2^2 + 1))^ ( (r^i, 0) 1 h; (r/i, 0) 1 /i||(Ai, ^I^) 1 0)i 
^ {-} 

X ((r/2,0) 112-AV2,0) l/2||(A2,/i2) 10)i((Ai,/ii) 10;(/i2,A2) 10||(A,/i) 10)^ 



(2.2) ^P (^' p) ^=1 -^=0 ■5=0 J=0 
-^(Ai,m)(M2,A2)5i=0S2=0 



(12) 



where { — } = {/i, Z2, Ai,/ii, A2,/i2,P, A,/i}. (The summation over r/i and rj2 is dropped 
because the pseudo-SU(3) model is a QTioj theory and inter-shell couplings are not allowed.) 
Since the pairing interaction acts only between like nucleons, the proton (neutron) pairing 
operator Hp {Hp ) can be obtained by replacing the quantum numbers in the generic 
expression, Eq. ([T^) , by their proton (neutron) counterparts. And as already noted, since 
the proton (neutron) pairing interaction is a SU(3) scalar operator in the neutron (proton) 
subspace, (A^, /iy) = (0,0) ((A^, ji-,,) = (0,0)), the expression for its matrix elements 



simplify considerably as compared to those of a general proton-neutron interaction, see 
Eg. (p!OD ■ Results for the matrix elements of Hp and Hp in a proton-neutron basis are 
given in Appendix ^ 

4 Pairing and collective states 

Since the primary objective of this study is to assess the importance of pairing correlations 
in heavy nuclei as revealed within the context of pseudo-SU(3) model, it is of interest to 
study the influence of the pairing interaction on the wavefunctions themselves. Consider 
the normal-space Hamiltonian 

Hqp = I ^'Q" ■ "-Q^ + G^Hl + G.H^p + aK] + hJ^ (13) 

where ^Q°' denotes the real quadrupole operator, J the collective angular momentum 
operator, and K'j stands for a residual interaction that is designed to generate i^-band 
splitting (see p6|, |2^). Although the deformation driving ^Q"" ■ '"Q" term is not exactly 
equal to its pseudo-space image Q" ■ Q"^, the difference is small. Indeed, within the 
leading SU(3) irrep the induced corrections to energies and electromagnetic transition 



probabilities have been shown to be smaller than 1% |T^. Since the exact results for the 
normal-space — > pseudo-space transformation are known for Hp and Hp (see Eqs. 15, 
P!7D , and the K'j and J^ operators transform as scalars, 

Hqp = I g» ■ g'^ + G^H}, + G.H^p + aK] + hJ^ (14) 

for practical purposes. (If not explicitly indicated otherwise, the parameters have the 
numerical values x = 4.32 keV, a = 0.202 MeV, and b = 9.26 keV which are compatible 
with a realistic set taken from best-fit calculations for heavy deformed nuclei [p8| ). 

The results presented next refer to the relatively simple system of two protons and two 
neutrons in pseudo-shells 77^ = 3 and rj^, = 4, respectively. Systems with a larger number of 
particles have also been studied and were found to display the same qualitative behaviour 



as the 27r2z/ system. The results for these more complicated cases will be published in [^ 
which focuses on complementary characteristics of the pairing interaction. An advantage 
of the present study is that concerns associated with truncating the model space can be 
avoided altogther because the 27r2z/ system, though rich in structure, is small enough to 
be solved without invoking truncation measures. 

Before taking a closer look at the quadrupole <-> pairing strength relation, it is instruc- 
tive to consider the spectrum for the case of pure pairing, that is, when x = a = 6 = 
and Gt^ = G,^ = 1 MeV in Hqp. The results shown in Fig. |], (A), provide a check on 
the calculations since in this case analytic results for the excitation energies can be given, 
and, (B), illustrate the structure of the pairing spectrum for a proton- neutron system. 



Figure ^ 
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Analytical results for the pairing spectrum can be derived in a quasi-spin formalism [|12|, |I^ 
with identical particle states classified according to their seniority s, which is simply the 
number of particles not coupled to J=0 pairs. For the combined proton-neutron system, 
the excitation energy Eg^g^ is given by the sum 

G 



Es^s^ = --r E isa[sa-2{n^ + l)] + 2N^{n^ + l)-N:, 



TT, U 



where s^ denotes the seniority, A'^^- the number of nucleons, and fig- the number of levels 
in the rif, shell (tt -^ proton, v — > neutron). (Recall that fi = (77 + l)(r7 + 2)/2) for 
the 77-th major (normal or pseudo) shell.) Taking A^ = N^ = 2 with ?]^ = 3 and 
rjy = A one obtains a pairing gap of 10 MeV (15 MeV) for the energy required to break 
a proton (neutron) pair. If one renormalizes the spectrum with respect to the ground 
state one obtains the results shown in Fig. |^. The states which are ordered according 
to increasing angular momentum group into degenerate sets classified by the seniority 
quantum numbers (s^, Sv)- The dimension of each set, which is given on the top of each 
horizontal bar, indicates that only a few states contribute to the low energy spectrum. The 
result of combining this pure pairing seniority level scheme with a quadrupole dominated 
rotational spectrum to yield more realistic results is considered next. 

To investigate the role pairing plays in a description of collective phenomena, consider 
the A^^ = Ny = 2 system introduced above but now with the pairing and quadrupole 
interaction strengths both taken to be non-zero. Specifically, the quadrupole strength x 
was fixed at a realistic value and the pairing strength was varied from to 0.30 MeV. 
(For simplicity, the proton and neutron pairing strengths were set equal to one another, 
Gtt = G,y = G.) An estimate for realistic Gj^ and Gi, values can be obtained from various 
phenomenological formulas, for example, in [2^ one finds the result Gjr = ^ MeV and 



Gi, = ^ MeV. For our example of a rare earth nucleus with 2 protons (neutrons) in the 
pseudo shell ?]^ = 3 [rji, = 4), which derives from a normal shell with principal quantum 
number r]^ = A {t]^ = 5), one has A = 136 and hence G^ = 0.125 MeV {G^ = 0.169 MeV). 
(Although this system happens to have the proton and neutron numbers of 52^Xe84, it is 
important to remember that the purpose of this investigation is to study the influence 
of the pairing interaction on excitation spectra and not to model a particular physical 
system. This exercise should be considered a forerunner for forthcoming attempts at 
describing and predicting experimental data. Best fit calculations for 52^Xe84 as well as 



other rare-earth nuclei will be published elsewhere [p8 |. 
Figure |, ||, |5| 



Figures ^ ^, and ^ indicate absolute values of the amplitudes (upper plots) and intensities 
(lower plots) of the calculated Oi, 4i, and 81 eigenstates. These values (vertical axis) are 
plotted as a function of the basis state number (running from left to right) and the pairing 
strength (increasing from front to rear). The basis states have been sorted as a function of 
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the eigenvalue of the second order SU(3) Casimir invariant, C2 = A^ + A/i + /i^ + 3(A + 11) ■ 
Hence, in each figure the most deformed SU(3) irrep is associated with the index 1, which 
is on the far left, and the least deformed with the highest index on the far right. The 
pairing strength is indicated on the right, where for simplicity the strengths were set equal 
to one another, Gt, = Gy = G. 

For G = the quadrupole interaction dominates and the ground state wavefunction is 
comprised of the most deformed SU(3) irrep only; that is, no mixing occurs, and this result 
is angular momentum independent. As the pairing strength is increased, contributions 
from irreps that correspond to less deformed configurations grow, but not all basis states 
contribute in the same way to each eigenstate. Instead, the results suggest that there may 
be a distinguishable pattern that applies for each value of the total angular momentum. 
If this turns out to be the case and can be predicted apriori, it could lead to a prescription 
for an angular momentum dependent model space truncation. This subject is currently 
under investigation. Regardless of whether this is true or not, the point to note is that 
for realisitic values of the pairing strength (0.1 MeV < G < 0.2 MeV) there is significant 
basis state mixing for all values of the total angular momenta — the pairing interaction 
definitely breaks the SU(3) symmetry. 

Notice that the contribution of the leading irrep increases with increasing angular 
momentum L. Two complementary factors contribute to this effect: First of all, within 
any U{Vi) representation the number of SU(3) irreps that contribute to a particular L 
value decreases as L increases. Indeed, if (A,/i) is the leading SU(3) irrep then the 
value L = (A + /i) is unique and the probability of it occuring in the eigenstate with 
J = L = (A+/i) is unity (100%), independent of the nature of the interaction. Secondly, as 
the value of L increases, the probability of finding pairs that couple to angular momentum 
zero decreases because the formation of a zero-coupled pair subtracts angular momentum 
from the system. Another way of saying the same thing is that the formation of a pair 
effectively reduces the system to one with two less particles and unless one is near the 
mid-shell region, the leading SU(3) irrep of such a system has a lower maximum L value. 

The influence of the pairing strength G on the rotational ground state band of the 
Ntt = Niy = 2 example is shown in Fig. ^ 



Figure ^ 



With the same parameters as used in the above analyses, one clearly finds that with G 
increasing from O.OMeV to 0.3MeV the system moves from a pure rotational structure 
towards one that is like that of the pure pairing Hamiltonian, see Fig. ^ for comparison. 
For realistic values of G (0.1 MeV < G < 0.2 MeV) there is clearly competition between 
the quadrupole dominated rotational picture and the pairing dominated seniority scheme. 
Effective moments of inertia 6 for each of the excited levels in Fig. M are shown in 



Fig. J^ The values plotted were extracted from the spectrum using the simple formula 
9 = -2E- '^^^ f^^^ ^^^^ ^^^ moments of inertia decrease with increasing pairing strength 
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is to be expected because typical 2+ states of rotational bands lie at a few tens of keV 
while those of pairing- dominated vibrator configurations are usually in the MeV range. 
Note that for a fixed value of the pairing strength the magnitude of the decrease in the 
moments of inertia is less the higher the L value. Again, this is the expected result 
because states with high angular momentum as compared to those with low values resist 
the formation of pairs because each pair subtracts angular momentum from the system, 
thereby making rotation more difficult to generate. So while from Fig. ^ an increase in 
the pairing appears to affect states with high L more strongly than those with lower L 
values, the opposite holds. The former is an illusion created by the L{L + 1) multiplier 
that enters into the energy equation; the pairing interaction is most effective in states 
with low L values. This picture is consistent, for example, with the interpretation of 
backbending as being associated with the break-up of a pair. 

5 Summary and outlook 

The formalism that is required for carrying out shell- model calculations within the frame- 
work of the pseudo-SU(3) model, extended to include exphcitly both the spin and proton- 
neutron degrees of freedom, has been introduced. In particular, the pairing interaction, 
which couples different irreps of SU(3), has been expressed in terms of pseudo-space unit 
tensor operators, and general expressions for calculating the matrix elements of these 
tensors, and hence the pairing interaction have been given. (Another study will give a 
more detailed account of the pairing force and explore additional features that reflect on 
the importance of pairing correlations in heavy nuclei [|1^.) The presented results make 
liberal use of an extended technology for calculating SU(3) coupling and recoupling co- 
efficients |TB| and a recently released code for evaluating SU(3) reduced matrix elements 
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The influence of pairing correlations on the collective (rotational) properties of the 
A^^ = 2 (r/^ = 3) and Nj^ = 2 [t],^ = 4) system was then considered. The Hamiltonian 
that was used includes a. Q"^ ■ Q"" term, the pairing interaction, and residual J^ and Kj 
terms that serve to generate the J{J + 1) and if-band splitting. (The splitting generated 
by J^ is trivial because rotational invariance insures that J is a good quantum number, 
but since Kj is not an exact symmetry the splitting generated by Kj is non- trivial, and 
can compete, for example, with that generated by the pairing interaction.) For vanishing 
pairing strengths (G^r = Gjy = 0), this Hamiltonian has non- vanishing matrix elements 
between states of the same irrep only. In particular, for J = states neither the J^ and Kj 
terms contribute so only the Q" ■ Q"^ term enters and as a consequence the leading SU(3) 
irrep, which is unique, is the only contributor to the ground state. Upon increasing the 
pairing strength, however, this perfect SU(3) symmetry is broken and other irreps enter, 
especially in the ground state. Similar results are observed for the states of higher angular 
momentum, but to a lesser degree because for these the role of the pairing interaction is 
curtailed and the Kj term enters. 

This paper is the first in a series of contributions which are aimed at achieving a proper 



13 



quantitative description of heavy nuclei within the framework of the pseudo-SU(3) model, 
with pairing playing its proper role. In keeping with this plan, the general formalism for 
calculating the matrix elements of any interaction was introduced, with explicit results 
given for the pairing interaction. To gain a deeper understanding of the importance 
of short-range correlations in strongly deformed systems, it will certainly be useful to 
investigate the effect pairing has on other collective nuclear properties such as qr factors, 
B(M1) values, Kj mixing, B(E2) values, etc. Such studies will serve as the basis for a 
quantitative description of heavy nuclei within the framework of the pseudo-SU(3) model. 
Beyond these studies, however, lies an even bigger challenge, namely, the development of 
a comprehensive program aimed at the unification of shell-model and collective-model 
theories. It is our goal to be able to give a proper account of macroscopic phenomena in 
terms of a user-friendly and usable microscopic theory. The simplifications provided by 
the pseudo-spin concept are, of course, key to achieving this lofty objective. 
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6 Appendix 

A Proton and neutron pairing matrix elements 

Applying the general result given in Section 2 for SU(3) matrix elements to the proton 
pairing operator, see Eq.(|12]), yields: 

{{K[K]< (a;, /x;) , K[f:^]ai (k, ^'^) > p\\\ /.o k'v {5;, s^ s'- f 



\^p\ 



{A'',r[/7r]a7r (A^, /i^r) , N^[fy\au {Xy, jJ-iy) } p (A, /i) kL {Sj,, Su} S ; J) 



X Y^ {{X, p) kL,{Xo, po) 1^0= lLo = 0\\{X', p') k'L')p 

Ao,Mo,p 



E 

Ptt 



(A^, /i^) (Ao, fio) (A^, p'^) Pn 

{X,,p,) (0,0) {XI, Pi) 1 

(A, ;U) (Ao, Po) (A', ;U') p 

p 1 p' 

{-} 

(15) 

where { — } = {Ai, pi, A2, P2, Pu}- The triple barred reduced matrix element that enters 
here is defined in Appendix p. The expansion coefficients Cf^^ "^'^"^ ^^-^ are given by: 



ff P12 (\o,Po) 



= E?(-l)'^^'^^'((2^i + l)(2/2 + l))^ 

X ((r/^, 0)1/1; (r^^, 0)l/i||(Ai, /xi)lLi = 0)] 



X ((r7.,0) l/2;(r?.,0) 1/2||(A2,P2) li^2 = 0), 



X ((Ai,/ii) 10; (/i2,A2) 10||(A,/i) 10) 



P12 



(16) 



where the constraint rj^,-^ = rj^,^ = rj^, has been invoked because within a pseudo-SU(3) 
model framework the action of the pairing operator is restricted to one major shell. (See 



Chapter 7 of |]T5[ for a more complete story, including, among other things, numerical 
results for the expansion coefficients for the 77 = 3 case.) 
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A complementary result (vr <-> i/, etc.) can be used to calculate matrix elements of the 
neutron pairing operator: 

{{KiKW. (a;, /^;) , KifiK (K. f^l) } p'(A', /i') k'l' {5:, s'j s'; J' 



\H'^\ 



{A^7r[/7r]tt7r (K, /^tt) , Klfu^u {K, /iz.) } P (A, yU) kL {S^, S^} S ] J) 



2J+1 



X ^ {{X, p) kL , (Ao, fio) Ko = lLo = 0\\ (A', fi') k'L' )p 

Ao,A«0,P 



E 

Ptt 



(A.,/iO (0,0) (a;,^;) 1 

(Ajy, /ii/) (Ao, ^o) (A^, /ij^) Piy 

(A, /i) (Ao, jio) (A', /i') p 
P 1 p' 

1 Pl2 (Xo , ^J.o) / AT' I f']^' f \' ,,'\Q' \\\n-pl2 {Xotl-lo) Svo=0 



X E ^&T)T;iA.)(^^[/^]«:'(^- p'j^: iii^cr^.o'jXA.f" III ^^^[^]«^(A., p.)s.), 

{-} 



(17) 



yP12 (Ao,A«o) 



with { — } = {Ai, /ii, A2, P2, Pu}- The C^(Ai u°)Tu2 A2) ^^ ^^^^ '^^^^ ^^ ^°^ neutrons, that is, 
r]^ -^ 7]y in Eq.([T|). 
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B Tripple-barred reduced matrix elements 

The Wigner-Eckart theorem for SU(2) gives 

( (A3, /X3) /ts/smsl T("^'^=^)'^^'^"^^ I (Ai, /ii) K.hmi ) 

= ( hmi, l2m2 I km, ) ( (A3, /ia) /^s^sl I T^""'' ^'^ ^^'^ \ \ {X,, fi,) K,h , ) (18) 

where ( (A3, ji^) f^shW T^'^^'^'^^'^^^^ \\ (Ai, /ii) /ti/i ) is a double-barred SU(2)-reduced matrix 
element. Analogously, it is possible to extend this factorization to SU(3) by introduc- 
ing SU(3)-reduced triple-barred matrix elements. An important difference in this case, 
however, is a sum over the multiplicity index p: 

( (A3, /i3) '^skms I T'^^^'^^)-^^^-^^ I (Ai, /iO Kikrm ) 

= J2 ((Ai,/ii)Ki/imi; (A2,/i2)/t2/2"^2| (A/i) /t3/3m3)p ( (A3, /i3) 1 1 1 t(-^2,/^2) III (;s^^^ ^^)^^ 
p 

= ^ {hmi, 121712 \hm3){(\i, fii) Kih; {\2, 1^2) ii2h\\ {hl^-i) ^.h) p 
p 

x((A3,/i3)|||T(^^''^^Mll(Ai, /ii))p. 
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Figure captions 



Figure 1: The group chains that are used for the classification of pseudo-SU(3) wave- 
functions are depicted in the two center columns. The eigenvalues of the corresponding 
Casimir invariants, which are indicated on the far left and right of the diagram, provide 
a complete labeling scheme for pseudo-SU(3) basis states. 



Figure 2: The excitation spectrum of Hqp for the pseudo- space configuration 
[{fp)^''^'^{gds)^"^'^] for the case of a pure pairing interaction, x = ^"^^ G^, = G^ = 1 
MeV. The horizontal axis denotes the total angular momentum with L = J since S* = 0. 
The seniority quantum numbers, {s^^, s,y) for protons and neutrons, respectively, are indi- 
cated on the far right and the numbers on top of each level bar denote the corresponding 
dimension. 



Figure 3: Spreading of the L = ground state over the 29 dimensional L = pseudo- 
SU(3) model basis of the [{fp)'^''^'^{gds)'^'^^'^] system. Absolute values of amplitudes 
(upper graphs) and intensities (lower graphs) are plotted on the vertical axis as a function 
of the basis state number (sorted according to the eigenvalue of the second order Casimir 
invariant of SU(3), C2 = A^ + A/i + /i^ + 3(A + ;u), and therefore according to their intrinsic 
deformation, with values decreasing from left to right) along the horizontal axis. The 
strength G of the pairing interaction is given on the far right, increasing front to rear 
from 0.0 to 0.3 MeV. For simplicity the proton and neutron pairing strengths were set 
equal, G = Gt, = G^. 
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Figure 4: Spreading of the L = 4 member of the ground band over the 101 dimensional 
L = 4 pseudo-SU(3) model basis of the [ (/p)^''^^(5'c?s)^''"^] system, see Fig. 3. 



Figure 5: Spreading of the L = 8 member of the ground band over the 101 dimensional 
L = 8 pseudo-SU(3) model basis of the [{fp)^''^^{gds)'^''^'^] system, see Fig. 3. 



Figure 6: Influence of the pairing strength (horizontal axis, G = Gj^ = G^) on members 
of the ground state of the [ {f p)'^^^'^ (gds)'^"^'^ ] system. Note that the rotational band on 
the left (G = 0) converts into a pairing dominated level structure (seniority level scheme) 
for {G = O.SMeV) on the right. A realistic range for the pairing strength is approximately 
(O.lMeV < G < 0.2 MeV). 



Figure 7: Moments of inertia as a function of the pairing strength {G = G^^ = Gy) for 
yrast states of the [{fp)^"'^'^{gds)^''^'^] system. The moments of inertia 9 were extracted 
from the calculated spectrum, see Fig. 6, by using the schematic formula 6 = ^^^ ' . 
Values of the angular momenta {J = L since S = 0) are indicated on the right. 
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